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Abstract 



The magnetic interference pattern and the spontaneous flux in uncon- 
ventional Josephson junctions of superconductors with d + is symmetry are 
calculated for different reduced junction lengths and the relative factor of the 
d and s wave components. This is a time reversal broken symmetry state. 
We study the stability of the fractional vortex and antivortex which are spon- 
taneously formed and examine their evolution as we change the length and 
the relative factor of d and s wave components. The asymmetry in the field 
modulated diffraction pattern exists for lengths as long as L = lOAj. 
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I. INTRODUCTION 



In the past several years one of the main questions in the research activity on high-T c su- 
perconductors has been the identification of the order parameter symmetry Jj]-|5[]. The most 
possible scenario is that the pairing state is an admixture of a dominant <i-wave with some 
small s-wave component. This fact is a direct consequence of the orthorhombic distortion 
of the systems which makes both the <i-wave and s-wave indistinguishable (they transform 
according to the identity representation of the group). There is a basic difference in the 
physics if one takes into account the phase difference between the two parts of the order 
parameter. The mixing due to orthorhombicity predicts a d+s or equivalently d—s order pa- 
rameter. This has been analyzed within the Ginzburg-Landau framework, valid close to T c 
H . Experimental observation of this possibility has been clearly realized in photoemission 
experiments ]7) and the c-axis tunneling ||. 



In addition to the above work, calculations based in BCS weak-coupling theory |5,ID 
predict that a mixed symmetry is realized in a certain range of interaction. This state has 
the time-reversal symmetry T broken. This symmetry is realized in bulk calculations only 
as a consequence of the absence of any orthorhombic distortion (the Fermi surface is either 
circular or tetragonal in the particular examples) which favors a phase difference of 71/ 2 
between the two components as opposed to n in the presence of it. 

The situation becomes more complicated if we consider surface effects. The observation 
of fractional vortices on the grain boundary in YBa2Cu307 by Kirtley et al. [TT|, may 
indicate a possible violation of the time-reversal symmetry near grain boundary (because 
the boundary breaks the bulk orthorhombic symmetry). Therefore it is interesting to study 
more this symmetry in the case of interfaces. 

In the present paper we study the static properties of one dimensional junction which 
contains a twin boundary where the pair transfer integral between the two superconductors 
has an extra relative phase in each twin. The maximum current I c that a junction can 
carry versus the external magnetic field H in direction parallel to the plane of the junction 
is calculated by solving numerically the Sine-Gordon equation. The stability of fractional 
vortices or antivortices f av which are spontaneously formed as a consequence of the 
symmetry, is examined in the absence of current and magnetic field for different lengths and 
relative phases. 

In the T-violated state the magnetic interference pattern as has been obtained by Zhu et 



al. 1 15] in the short junction limit is assymetric. They conclude that for a long junction the 



magnetically modulated critical current is basically identical to the conventional 0-0 junction 
due to the formation of the spontaneous vortex near the center of the junction. Our exact 
numerical calculations, show that there is a "dip" near the center of the diffraction patterns 
even for junctions as long as 10Aj. 

The rest of the paper is organized as following. In section II we discuss the Josephson 
effect for a mixed wave symmetry. In section III we present the results for the magnetic flux 
and the interference pattern. Finally, summary and discussions are presented in the last 
section. 
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II. JOSEPHSON EFFECT FOR A MIXED WAVE SYMMETRY 



We discuss the Josephson coupling at the interface between two superconductors (A and 
B) both with a two component order parameter (nf , nf )■ We can think of the interface 



as a Josephson junction, so the Josephson current phase relation is |L2| 



J= £ J^sin(0f-0/), (1) 
«,i=i 

where J C ij is the coupling between the components nf on side B and with nf on side A 
[nf = \nf \exp(i(pj)]. We consider some special cases. 

(i) For (i-wave symmetry one component of the order parameter vanishes at the interface 
(n 2 = 0). The Josephson current density becomes J = | J c ii| sin(0 + ir) with J cll < 0. 

(ii) For d + s-wave we are restricted to the case where <f>f — <fif = (f)f — (j)2 = it is fixed 
on both sides of the interface. The current density J depends only on one phase difference 
through the interface, say (ft = (fif — (pf 

J(0) = |J c |sin(0 + #) (2) 



Jc — Jell + Jc22 — Jc\2 — Jc21 (3) 

with 9 = for J c > and 9 = n for J c < 0. 

(Hi) For d + is-wave case the intrinsic phase difference within each superconductor A 
and B can be assumed to be (fif — 0^ = 0f — 0f = ir/2. The current density J is 

J(0) = J c sin(0 + #) (4) 

with 



= vVdl + ^c2 2 ) 2 + (Jcl2 - ^c2l) 2 , (5) 

tan(^) = - — — — ■ (6) 

^cll + ^c22 

We consider two superconducting sheets A and C which overlap for a distance L with the 
superconducting sheet B, in the x-direction, (as shown in Fig 1). All three superconducting 
sheets have dominant (i-wave symmetry, with a small s-wave component. The angles between 
the crystalline a axis of each superconductor A, B, C with the junction interface are defined 
as #i,#2,#3- We describe the entire junction with width w small compared to Aj in the y 
direction, of length L in the x direction, in external magnetic field H in the y direction. 
The intrinsic phase difference 9(x) is <p c i in < x < |f and <p C 2 in |f < x < L. The 
superconducting phase difference <fi across the junction is then the solution of the Sine- 
Gordon equation 

-^ = —8m[<i>(x) + e(x)], (7) 
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with the inline boundary condition 



~T \x=0,L = ±7T + H 

dx 2 



(8) 



The Josephson penetration depth is given by 



A., 



he 2 



\ 8nedJ c 



where d is the sum of the penetration depths in two superconductors plus the thickness of 
the insulator layer. We also assume that J c is constant within each segment of the interface. 

To check the stability we consider small perturbations u(x,t) = v(x)e st on the static 
solution <f>(x), and linearize the time-dependent Sine-Gordon equation to obtain: 



d 2 v(x) 
dx 2 



+ cos[0(x) + 9(x)]v(x) = Xv(x), 



(9) 



under the boundary conditions -^-\ x =o,l = 0, where A = — s 2 . It is seen that if the 
eigenvalue equation has a negative eigenvalue the static solution <fi(x) is unstable. 

We can also compute the free energy of the solution for zero current and external magnetic 
field 



TiJ c w r L 
2e Jo 



cos [0(x) + 6{x)} + ^ (?£ 



(10) 



Note that the no vortex solution = everywhere is not a solution of this problem. 

When (p c i = (f> C 2 = we have the conventional s-wave junction. In case <p c i = 0, 4> C 2 = ir 
we have the c?-wave or d + s-wave junction. The above cases have time reversal symmetry 
(T-conservation). When cl , c2 are slightly different from and 7T, we have the d + is- wave 
pairing, which is a broken time reversal symmetry state (T- violation). In this work, the 
particular parameters we use are c i = 0.017T, <p C 2 = 1.087T, and the pairing state is d + is. 



III. SPONTANEOUS MAGNETIC FLUX AND INTERFERENCE PATTERNS 
FOR THE T- VIOLATING PAIRING STATE 

In Fig. 2 we plot the maximum current for a symmetric — it junction as a function 
of the magnetic flux $ (in units of $o — ff ) for different junction lengths: (a) L = 10, (b) 
L = 4, (c) L — 2, (d) L = 1 (A j = 1). The circles and squares in this figure correspond to 
the fractional vortex (/„) and antivortex (f av ) branch. For most of the range of existence of 
fv (fav) the magnetic flux is positive (negative) while there is a small region where it turns 
into antivortex (vortex). Similar calculation has been done ||13|| , |l4jl who considered the f v , 



since in this case the plot is symmetric in H. As we can see, there is a "dip" at $ = 0, 
for lengths as long as L = 10. In Fig. 3 we present our calculations for the T- violation 
case where <p c \ = 0.017T and = 1.087T. We also plot for L = 1 the analytical result (solid 



line) of Zhu et al. [15]. In contrast to the pure c?-wave case, for small lengths this pattern is 
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assymetric and the "dip" in the maximum current does not occur at $ = 0, but at a finite 
$ value. This behavior also exists for lengths as long as L = 10. 

If we plot I c vs H (and not $) then the two branches in Fig. 3, will be almost coinci- 
dent and one might draw the conclusion that the behavior for a long junction is the same 
independent of the symmetry. The proper quantity to consider though is the total magnetic 
flux which includes both the contribution from the external field and the induced self field. 
It should be remarked that for an s-wave junction the relation between $ and H is linear 
for small H so that the plot of I c vs if or $ does not show any differences for small H. 
For higher H however the overlapping branches (for long L) are unfolded. In the case of a 
different symmetry even the small H form can change due to the existence of spontaneous 
magnetization. In this case, if we consider the zero current solutions, and vary the magnetic 
field, both f v , f av are stable, whereas for L = 4, L = 2, L = 1 the stable regions in the 
magnetic field are separated by the unstable ones, and this behavior persists in both the /„, 
f av cases. Also in the short junction limit (i.e. L = 1) these two branches coincide at the 
maximum current. 

Figure 4 addresses the question of spontaneous flux generation in junctions with broken 
time reversal symmetry (T-violation) as a function of the reduced length (L) and the relative 
factor of s and d components. The long dashed line is the result of [|T5[] which compares with 
our numerical result (solid line). Both cases have <p c i = O.Oi7r,0 c2 = 1.087T. We have also 
used two other values for c2 i.e. 0.97T (dotted line) and 0.87T (dashed line). We conclude 
that as we decrease the value of c2 the fractional vortex /„ tends out to be a 27r vortex, 
whereas the fractional antivortex gradually loses its flux content. 

In Fig. 5a we have plotted the magnetic flux $ (solid line, $0 — 1) versus the value c2 
for L = 10 and H = 0. We see that as we decrease c2 the magnetic flux increases linearly 
for the f v branch from <3> « 0.45, for c2 = 1.087T to $ ~ 1 for c2 = 0. But this last point is 
unstable, as seen by the stability analysis from which the lowest eigenvalue is also displayed 
(light line). This is expected since it goes to a point in the unstable (1,2) branch of the 
usual s-type junction fi~6fl . The /„„ branch increases linearly its flux as we decrease c2 and 
goes to the stable (0, 1) branch. Here we follow the notation of |17[]. This linear dependence 
of $ from (f) c2 can also be seen in the analytical result of [[15] for large lengths, where the 
approximation they made is valid. On the other hand as we increase </> c2 from 1.08-7T to 2ir, 
the f v branch decreases its flux and goes to the stable (0, 1) branch, while the f av branch 
goes to the unstable (—2, —1) branch. When we change <f> c \ and keep c2 = 0, from to 
2tt, the (0, 1) branch goes to f v and then to the unstable (1, 2), while the unstable (—2, —1) 
branch goes to f av and then to the stable (0, 1). The situation is a little bit different for small 
lengths as can be seen from Fig. 5b where L = 1, H = 0. Here f av branch is unstable for 
4> C 2 = 1.08-7T and by decreasing c2 it gets stabilized, but f v branch is stable for c2 = 1.087T 
and then it becomes unstable. We can see that from Fig. 6 where we plot the ratio F j Fq of 
the free energy of the state with some spontaneous flux to the state with no flux. This ratio 
becomes larger than one as we decrease the c2 , for the /„ branch, for small lengths. On the 
other hand, when F/F < 1 the no flux state is metastable and the final state will be the 
one with spontaneous flux. Notice that the magnetic flux remains almost constant - almost 
zero - which can be expected since we are in the short junction limit where self currents are 
neglected. 
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IV. CONCLUSIONS 



We have studied the static properties of a one dimensional junction with d + is order 
parameter symmetry. The magnetic interference pattern is asymmetric, and there exist a 
"dip" near $ = for lengths as long as lOAj. The diffraction pattern of a junction can give 
us information about the pairing symmetry, at least where junctions are formed. 

We have followed the evolution of spontaneously formed vortex and antivortex solutions 
for different mixing between the s and d components of the order parameter. We have shown 
that for small lengths the fractional vortex becomes unstable as we decrease the extra phase 
of the pair transfer integral in the right part of the junction. We conclude that when a mixing 
state symmetry is realized, the fractional vortex and antivortex solutions evolve differently 
and this characterizes the d + is-wave pairing. We expect these findings to hold even if a 
bulk d + s state evolves continuously as a function of distance from the interface to a d + is 
one, as long as the there is a well defined area close to the interface where the time reversal 
symmetry is not conserved and the junction is formed. 
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FIGURES 



FIG. 1. The one dimensional junction geometry. The dashed line marks the twin boundary. 

FIG. 2. Critical current I c versus the magnetic flux $ (in units of 3>o) for a symmetric — ir 
junction, for different junction lengths: (a) L = 10, (b)L = 4, (c)L = 2, (d)L = 1. 

FIG. 3. Critical current I c versus the magnetic flux $ for a junction with d + is symmetry, 
4> c i = 0.017T, 4> c2 = 1.087T, for different junction lengths: (a) L = 10, (b)L = 4, (c)L = 2, (d)L = 1. 

FIG. 4. The spontaneous magnetic flux $ as a function of the reduced junction length L, for 
different values of the intrinsic phase (p C 2 in the right part L/2 < x < L of the junction and 
0ci = O.Olvr. 

FIG. 5. The evolution of the fractional vortex f v and antivortex f av as a function of (p C 2, for 
two different lengths (a) L = 10, (b) L = 1. The stability of theses branches is also denoted by the 
lowest eigenvalue Ai of the linearized eigenvalue problem. Note that the double arrow connects the 
flux with its stability curve. 

FIG. 6. The ratio of the free energy F/Fq as a function of the reduced junction length L, for 
different values of <p C 2- (a) fv, (b) fav 



8 



x=0 x=L/2 x=L 



T - 
1 


► 


A | 

1 


C 




















B 




I 


— ► 



x 



2.0 




o.o 

2 4 6 8 10 

L 



